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INVESTIGATION  OE  THE  WORKING  CASCADES  OP  TURBINES  AT 
SUPERSONIC  VELOCITIES 


A.  V.  Gubarev 


The  possibility  of  developing,  in  a  single  stage  of  a  turbine, 
large  heat  drops  without  a  noticeable  decrease  In  the  efficiency  Is  a 
very  attractive  Idea.  However,  up  to  now  there  have  been  no  reliable 
methods  for  profiling  working  cascades  which  operate  with  satisfactory 
efficiency  at  high  supersonic  velocities.  In  addition,  up  to  now 
there  have  been  no  Intense  studies  of  working  cascades  at  such  veloci¬ 
ties. 

In  this  paper  we  present  certain  results  of  an  investigation  of 
working  cascades  at  supersonic  velocities  [l].  Profiles  of  the  work¬ 
ing  cascades,  built  according  to  the  recommended  methods,  made  it 
possible  to  noticeably  Increase  the  efficiency  of  certain  turbines  of 
the  Kaluga  factory. 

1.  Features  of  Supersonic -Flow  Structure  In  Working 
Cascades  and  the  Method  of  Designing  New  Cascades 

In  the  literature  [2]  it  is  usually  noted  that  at  supersonic 
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velocities,  thick  leading  edges  of  the  airfoils  of  a  cascade  are  the 
source  of  additional  losses.  However,  recent  studies  [5,  4]  have 
shown  that  under  certain  conditions,  tapering  of  the  leading  edge  does 
not  produce  the  desired  result.  Moreover,  Oswatlsch  [3]  has  established 
that  even  In  the  case  of  a  knife-edge  on  the  airfoils  a  complex  system 
of  head  shocks  forms  at  the  cascade  Inlet.  An  analysis  of  the  results 
obtained  at  the  Moscow  Power  Instltue  [2,  5]  and  those  obtained  by 
other  authors  shows  that  the  head  shocks  occurring  ahead  of  the  cascade 
can  be  divided  Into  three  groups: 

1.  those  caused  by  flow  past  the  leading  edge  —  "edge  head 
shocks, " 

2.  those  that  depend  on  the  shape  of  the  vane  channel  —  "choking 
of  the  vane  channel, "  and 

3.  those  caused  by  an  off-deslgn  flow-inlet  angle  —  an  "oblique 
frontal  shock. " 

Occurring. ahead  of  the  cascade,  they  combine  and  form  a  complex 
system  of  head  shocks;  therefore  they  cannot  be  examined  Individually 
during  experiments.  However,  such  separate  examination  Is  expedient, 
since  It  would  make  It  possible  to  find  a  way  to  Improve  working  air¬ 
foil  cascades  designed  for  high  supersonic  velocities. 

To  decrease  the  Intensity  of  the  "edge  head"  shocks  we  must 
attempt  to  decrease  the  thickness  and  the  design  angle  of  the  leading 
edge.  However,  such  measures  do  not  always  make  It  possible  to  decrease 
the  Intensity  of  the  head  shocks.  This  Is  observed  when  a  "choking" 
shock  forms  ahead  of  the  cascade  and  a  subsonic  stream  flows  past  the 
leading  edges.  We  should  mention  that  "edge  head"  shocks  are  local 
In  nature,  while  "choking"  shocks  always  partition  off  the  entire 
flow.  Therefore  we  can  state  that  a  great  percentage  of  the  losses 
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In  head  shocks  Is  due  to  the  presence  of  "choking"  shocks. 

Consequently,  to  decrease  the  losses  In  head,  shocks  we  must,  in 
addition  to  decreasing  the  thickness  and  design. angle  of  the  leading 
edge,  also  design  the  vane  channel  of  the  cascade  such  that  "choking" 
shocks  do  not  occur  [6]. 

To  decrease  the  intensity  of  "ohllque  frontal  shocks"  the  "effec¬ 
tive"  slope  of  the  inlet  sector  of  the  hack  of  the  airfoil  must  he 
close  to  the  flow-inlet  angle  into  the  cascade  [5]*  We  should  note 
that  the  problem  of  lowering  the  intensity  of  the  head  shocks  at  the 
cascade  inlet  should  be  solved  complexly,  since  the  above  requirements 
are  mutually  exclusive  in  certain  cases. 

Let  us  examine  certain  features  of  supersonic  flow  in  the  vane 
channel  of  a  cascade.  We  will  consider  that  there  are  no  head  shocks 
and  that  the  flow  is  uniform  at  the  cascade  inlet. 

Since  disturbances  in  supersonic  flow  are  propagated  only  along 
the  stream,  the  flow  past  the  back  of  the  airfoil  in  the  inlet  section 
will  be  the  same  as  flow  past  an  isolated  body.  If  this  part  of  the 
back  of  the  airfoil  is  curvilinear,  as  is  the  case  in  cascades  for 
subsonic  velocities,  there  will  be  intense  flow  acceleration  along  it. 
Such  overexpansion  of  the  flow  results  in  the  formation  of  an  intense 
compression  shock  and  a  turbulent  zone  about  the  back  of  the  airfoil 
[ji,  1].  According  to  Oswatlsch  [3],  the  losses  in  the  turbulent  zone 
reach  6-8^.  Therefore  we  must  exclude  the  possibility  of  flow  accel¬ 
eration  in  the  inlet  section  of  the  airfoil.  This  problem  can  be 
solved  as  follows. 

First  method  —  the  airfoil  back  in  the  inlet  section  is  made 
with  reverse  concavity.  This  is  the  method  of  "stepwise  flow  stagna- 
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Second  method  —  flow  stagnation  In  a  steep  shock  occurring  at  the 
leading  edge  from  the  concave  surface. 

In  both  cases  the  vane  channel  Is  convergent-divergent,  or  diver¬ 
gent,  since  after  the  head  shocks  at  the  Inlet  section  of  the  vane 
channel  the  flow  is  stagnant. 

Supersonic  cascades  are  usually  used  In  stages  with  low  volume 
gas  passage;  therefore  they  should  be  designed  such  that  the  total 
losses  are  minimum.  As  has  been  shown  by  studies  at  the  Dzerzhlnskly 
All-Union  Heat  Engineering  Institute  and  the  Moscow  Power  Institute, 
for  this  there  must  be  relatively  low  flow  rates  at  sections  of  maximum 
curvature  In  the  vane  channel.  Consequently,  the  terminal  and  airfoil 
losses  In  cascades  operating  at  supersonic  velocities  can  be  decreased 
by  the  same  method  —  stagnation  of  the  flow  in  the  Initial  section  of 
the  vane  channel.  However,  the  degree  of  stagnation  of  the  flow,  from 
the  standpoint  of  lowering  the  total  losses,  can  be  greater  than 
necessary  to  Insure  minimum  airfoil  losses. 

Figure  1  shows  types  of  airfoils  designed  according  to  recommenda¬ 
tions.  We  investigated  about  JO  different  cascades  that  had  differing 
airfoil  shapes,  relative  pitch,  and  angles  of  Incidence.  The  table 
gives  the  basic  geometric  characteristics  of  the  cascades;  we  will 
next  discuss  the  results  of  their  investigation. 

2.  Results  of  Investigation  of  the  Cascades 
A.  Investigation  of  Type  V  Cascades 

The  geometric  features  of  type  V  airfoils  (Fig.  la)  are  as 
follows;  the  thickness  and  design  angle  of  the  leading  and  trailing 
edges  are  small;  the  elongated  sections  on  the  airfoil  back,  at  the 
flow  Inlet  and  In  an  oblique  section,  and  small  sectors  on  the  convex 


surface,  near  the  edges,  are  straight;  the  curvature  of  the  hack  and 
the  convex  surface  near  the  vane  channel  changes  smoothly. 

Figure  2  gives  the  curves  of  pressure  dlstrlhutlon  along  the 
airfoil  TR-iV  in  cascade  No.  2.  The  vane  channel  in  this  cascade  is 
convergent-divergent,  with  its  mlnimim  cross  section  located  at  the 
beginning  of  the  channel  (points  7-8  and  15-14).  Examination  of  the 
curves  shows  that  with  supersonic  flow-inlet  velocities  (M  >  l),  a 
plane  shock  occurs  ahead  of  the  cascade  channel;  with  increasing  speed 
this  shock  shifts  along  the  stream.  For  regimes  M  >  1.  5  the  head 
shocks  enter  the  vane  channel.  In  it  the  flow  is  accelerated,  and  at 
subsonic  velocities  in  the  minimum  cross  section  its  deep  overexpansion 
is  noted  at  the  back  and  at  the  convex  surface.  At  supersonic  veloci¬ 
ties  the  minimum  dimensionless  pressure  is  decreased  considerably, 
and  shifts  noticeably  along  the  flow.  The  Intensity  of  the  diffuser 
section  at  points  14-17  decreases. 

Figure  2b  gives  a  comparison  of  the  pressure  distribution  along 
the  airfoil  in  cascade  No.  4  for  three  flow- inlet  angles.  The  vane 
channel  of  this  cascade  has  two  narrow  sections,  one  at  the  inlet 
and  the  other  at  the  exit  from  the  channel.  After  the  second  throat 
the  channel  expands.  An  examination  of  the  curves  shows  that  in  sub- 
sonl,  regimes  the  pressure  distribution  along  the  airfoil  depends 
little  on  the  flow-inlet  angle.  At  supersonic  velocities  (M  =  1.45) 
a  change  in  the  flow- inlet  angles  within  the  same  range  results  in  a 
substantial  change  in  the  pressure  curves:  with  increasing  angle,  the 
flow  velocity  at  the  inlet  section  of  the  back  of  the  airfoil  decreases, 
with  a  simultaneous  abrupt  increase  in  the  flow  stagnation  intensity 
in  the  compression  shock. 

Such  a  change  in  the  pressure  curves  can  be  explained  as  follows: 
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Vflth  an  Increase  In  the  flow-inlet  angle  ahead  of  the  cascade  there 
occurs  an  "oblique  frontal  shock,"  and  the  flow  velocity  at  the  Inlet 
decreases;  the  Increased  Intenaity  of  flow  stagnation  at  the  Inlet 
to  the  vane  channel  Is  caused,  in  turn,  by  the  displacement  of  the 
"choking"  shock  against  the  flow,  caused  by  decreased  velocity. 

Figure  3  gives  the  curves  of  the  dependence  of  the  total  losses 
In  the  cascades  of  airfoils  TR-IV  and  TR-2V.  In  the  case  of  convergent- 
divergent  vane  channels  (Plg.  3a)  the  maximum  losses  are  reached  at 
M  =  1.0.  The  losses  decrease  in  subsonic  and  supersonic  regimes. 
Characteristically,  maxlmiun  losses  at  M  =  1. 0  are  found  In  cascade 
No.  1  (at  f  =  0.657),  while  the  minlmtun  losses  occur  In  cascade  No.  5 
(at  f  =  0.909) •  This  fact,  and  also  the  nature  of  the  change  of 
losses  versus  M,  agree  with  data  from  Laval-nozzle  studies.  With 
Increased  flow-inlet  angle  the  losses  decrease  throughout  the  entire 
range  of  velocities,  which  can  be  explained  as  follows:  first,  for 
large  Inlet  angles  flow  stagnation  occurs  stepwise  (in  an  "oblique 
frontal  shock  and,  subsequently,  at  the  leading  edge  and  ahead  of  the 
vane  channel),  and  second,  in  the  investigated  types  of  cascades  the 
channel  at  the  Inlet  Is  divergent,  and  therefore  there  Is  observed  no 
Increase  In  the  flow-outlet  angle  for  high  M  values  of  flow  expansion 
in  the  oblique  section  of  the  cascade  [5]- 

Figure  3b  gives  curves  of  the  total  losses  In  cascades  of  airfoils 
with  a  slight  change  In  the  width  of  the  channels.  This  type  of  cas¬ 
cade  Is  characterized  by  weak  dependence  of  losses  on  Mach  number. 

Only  at  =  17°  (cascade  No.  4)  do  the  losses  sharply  Increase  In 
M  =  1. 0  regimes.  Such  behavior  of  the  curves  Is  caused  by  the  struc¬ 
ture  of  the  shocks  and  the  expansion  waves  at  the  cascade  Inlet. 
Actually,  when  =  ±7°  an  expansion  wave  forms  ahead  of  the  cascade. 


-6- 


’  since  p.  <  The  flow  accelerates.  Therefore,  at  low  supersonic 

velocities  the  intensity  of  the  "choking"  shock  Increases  conslderahly 
compared  with  the  case  and  then,  for  high  Mach  niimbers,  when 

the  head  shocks  enter  the  vane  channel,  the  losses  decrease  and  even 
become  somewhat  smaller  than  for  large  flowrlnlet  angles.  This  can  be 
explained  by  the  fact  that  when  =  ±7°  there  is  no  flow  deviation 
in  the  oblique  section. 

B.  Investigation  of  Cascades  Consisting  of 
Modified  TR-IV  Airfoils 

To  explain  the  influence  of  the  geometric  shape  of  the  vane  chan¬ 
nel  and  the  airfoil  on  the  characteristics  we  investigated  Cascades 
Nos.  6-11.  The  TR-IV-II  airfoil  differs  from  the  TR-lV  in  that  there 
is  a  break  in  the  leading  section  of  the  back  and  in  the  middle  cross 
sections  of  the  vane  channel  the  convex  surface  and  the  airfoil  back 
have  arcs  of  greater  radius.  There  is  no  break  in  the  back  of  the 
leading  section  of  airfoil  TR-IV-I. 

From  an  examination  of  the  spectra  and  also  measurements  of  the 
pressure  distribution  along  the  profile  we  can  draw  the  following 
•  conclusions:  a)  at  low  supersonic  velocities,  flow  stagnation  at  the 

cascade  inlet  is  caused  by  channel  blockage;  b)  with  increasing  veloc¬ 
ity  the  shock  shifts  along  the  flow,  and  only  at  M  ^  1.  5  is  the  shock 
at  the  break  in  the  airfoil  back  (Fig.  4-a),  l.e.,  the  effect  of  step¬ 
wise  flow  stagnation  begins  to  be  manifested  at  the  cascade  inlet. 

For  greater  clarity,  graduation  lines  are  plotted  on  the  airfoil  sur¬ 
faces;  as  a  result,  the  structure  of  the  characteristic  is  quite 
evident . 

For  cascades  Nos.  6  and  9  it  is  characteristic  that  the  minimum 
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section  of  the  vane  channel  is  at  the  inlet  section.  Main  expansion 
of  the  channel  occurs  at  the  point  of  maximum  curvature.  In  Pig.  4b 
(flow  spectrum  for  cascade  No.  ,9)  we  see  the  shocks  in  the  vane  channel. 
For  a  broad  set  of  regimes  (M  =  1.2-1. 9)  the  position  and  shape  of  the 
shocks  did  not  change.  This  shock  Is  caused  by  oVerexpanslon  of  flow 
with  flow  past  the  airfoil  back  and  by  intense  compression  of  the  flow 
along  the  convex  surface  [5]. 

Figure  3c  gives  curves  of  the  dependence  of  losses  in  cascades 
consisting  of  TR-IV-I,  TR-IV-II,  and  TR-lV-2  airfoils.  Since  the  vane 
channels  of  these  casccides  were  divergent  at  the  outlet,  we  detect  the 
characteristic  dependence  of  losses  on  Mach  number:  when  M  =  1  there 
are  maximum  losses,  which  then  decrease  considerably.  For  cascades 
Nos.  6  and  9  the  calculated  regimes  (M  =  1.6),  determined  from  the 
ratio  of  the  outlet  area  to  the  minlmiim  area,  were  not  reached  during 
the  experiment,  and  therefore  the  losses  continually  decreased  with 
increasing  velocity.  We  should  note  that  when  M  =  0. 5-I.3  the  losses 
are  greater  In  cascade  No.  S,  while  when  M  >  I.3  the  losses  are  higher 
In  cascade  No.  6  (by  2J^);  cascade  No.  6  differed  from  No.  9  only  In 
the  absence  of  a  break  in  the  outlet  section  of  the  airfoil  back. 
Consequently,  stepwise  stagnation  Is  expedient  at  high  supersonic 
velocities.  * 

The  high  level  of  losses  (about  17^)  In  cascades  Nos.  7  and  10  In 
the  calculation  regimes  (M  =  1.4-5)  can  be  explained  by  the  fact  that 
most  Intense  expansion  of  the  channel  occurs  before  Its  maximum  curva¬ 
ture.  Therefore  the  flow  at  the  airfoil  back  Is  overexpanded,  and  an 
Intense  shock  occurs  In  the  channel;  as  observations  showed.  Its 
intensity  was  higher  than  In  cascade  No.  9  (see  Pig.  4h).  In  addition, 
the  secondary  flows  in  this  cascade  were  also  relatively  great. 
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In  cascade  No.  11  (airfoil  TR-lV-2)  the  losses  at  M  =  0.5-IO 
were  lower  than  In  cascades  No.  6-10,  which  can  be  explained  by  the 
high  value  of  the  parameter  f  and  the  position  of  the  maximum  cross 
section  at  the  outlet  section  of  the  channel.  Minimum  losses  (about 
13.5^)  are  attained  at  regimes  M  =  1.35-1»55>  which  are  the  calculation 
regimes  for  the  given  cascade. 

C.  Investigation  of  VS-Type  Cascades 

The  geometric  features  of  type  VS  airfoils  (Plg.  lb)  are  as 
follows:  a)  a  convex  surface  formed  by  an  arc  of  a  single  radius;  b) 
the  design  angles  of  the  leading  and  trailing  edges  are  zero.  The 
TR-lVS-2  airfoil  has  a  break  (5  =  10°)  at  the  inlet  part  of  the  back. 

The  Inlet  sector  of  the  vane  channel  is  mlnlmm  and  the  channel  diverges 
uniformly  (f  =  O.827).  For  the  TR-IVS-3  airfoil 'the  break  Is  decreased 
to  5°  and  the  vane  channel  is  convergent-divergent  (7  =  0.9)*  The 
channel  width  Is  constant  in  a  cascade  consisting  of  TR-lVS-5  airfoils. 

Figure  4c  shows  the  flow  spectrum  of  cascade  No.  12.  At  low 
supersonic  velocities  we  detected  two  head  shocks:  a  disconnected  shock, 
ahead  of  the  break  point  on  the  back,  and  a  curvilinear  shock  ahead 
of  the  edge.  With  increasing  Mach  number  the  disconnected  shock  becomes 
oblique  while  the  curvilinear  shock  shifts  along  the  flow,  and  In 
regimes  greater  than  M  =  1.5  Its  right  branch  enters  the  vane  channel. 
Flow  stagnation  at  the  cascade  inlet  occurs  in  a  system  of  oblique 
shocks. 

Figure  3c  gives  curves  of  the  total  losses  In  VS  cascades.  The 
minimum  losses  In  these  cascades  correspond  to  the  calculation  values 
of  M  for  the  divergent  part  of  the  channel.  Characteristically,  when 
M  >  M  the  losses  in  cascade  No.  I3  Increase  only  Insignificantly. 
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Analogous  behavior  of  the  curves  was  also  observed  for  cascades  Nos. 

10  and  11.  The  losses  In  No.  14  for  M  >  1.3  are  higher  than  In  Nos. 

12  and  I3. 

D.  Investigation  of  the  Plow  Field  Behind  the  Cascades 

Figure  5  gives  certain  characteristic  curves  of  the  loss  distri¬ 
bution  along  the  height  of  the  cascades.  In  subsonic  regimes,  as  a 
x*ule,  there  are  Intense  loss  peaks  caused  by  the  secondary  flows.  With 
Increasing  M  the  loss  peaks  decrease;  however,  for  most  of  the  Investi¬ 
gated  cascades  the  losses  in  the  average  cross  section  (airfoil  losses) 
Increase.  This  Is  due  to  boundary  layer  separation  In  transonic  regimes. 
With  a  further  Increase  In  velocity  (M.  >  l)  most  cascades  are  character¬ 
ized  by  a  considerable  decrease  of  the  loss  peaks  near  the  end  walls 
and  also  airfoil  losses.  At  high  supersonic  regimes  the  end  losses  In 
the  cascades  are  very  slight. 

Figure  6  gives  distribution  curves  for  flow-outlet  angles  along 
the  height  of  the  cascade.  With  Increasing  M  nonuniformity  of  flow 
along  the  height  of  the  cascade  decreases  considerably,  which  Is 
reflected  In  a  decrease  In  the  Intensity  of  secondary  flows  In  the 
cascades.  We  must  also  note  that  In  all  the  Investigated  Impulse  cas¬ 
cades  there  was  no  flow  expansion  In  the  oblique  section  for  regimes 
M  >  1. 

Conclusions 

The  results  of  systematic  studies  conducted  at  the  Moscow  Power 
Institute  showed  the  expediency  of  designing  three  different  types  of 
working  cascades:  for  subsonic,  transonic,  and  supersonic  velocities 
[7,  8].  Prom  our  work  we  can  draw  the  following  conclusions: 

1.  We  have  verified  the  developed  recommendations  on  the  design 
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of  more  effective  cascades  for  supersonic  velocities  (loss  levels  1.5- 
2-times  .lower).  However,  the  loss  level  (13-15/^)  should  not  yet  be 
considered  satisfactory. 

2.  The  efficiency  of  working  cascades,  at  transonic  velocities 
and  the  calculation  Mach  nvunber  are  determined  mainly  by  the  parameter 

^  =  ^mln/^2- 


3.  In  regimes  M  >  Mp  the  losses  in  the  working  cascades  increase 
Insignificantly,  which  allows  us  to  recommend.  Just  as  for  Laval  noz¬ 
zles,  designing  cascades  from  a  calulation  in  order  to  use  them  in 
regimes  greater  than  those  calculated.  However,  the  physical  signifi¬ 
cance  of  these  recommendations  Idffers  from  analogous  ones  proposed 
for  guide  cascades,  since  flow  is  not  deflected  in  the  oblique  section 
of  working  cascades. 

To  decrease  losses  in  the  head  shocks  for  high  Mach  nvimbers 
it  is  expedient  to  use  stepwise  flow  stagnation. 

5.  For  high  supersonic  velocities  the  end  losses  in  cascades 
developed  by  the  Moscow  Power  Institute  are  very  slight  compared  with 
the  total  losses.  This  makes  it  possible  to  recommend  the  use  of 
these  cascades  in  the  first  control  stages  of  turbines. 
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straight  section 


Fig.  1.  Shapes  of  airfoils  (Moscow  Power  Institute) 
for  supersonic  velocities:  a)  V-type  airfoil, 
b)  VS-type  airfoil. 


TABLE 


Cascade  nueber 

. 

Type  of 
airfoil 

Cord 

b 

Relative  pitch, 
t 

S 

9 

9 

Effective 
flow  outlet 
angle 

^seff 

Ibitio  of 
BiniJDun 
cross  sec¬ 
tion  to 
outlet  oros 
section 

•y _  “^in 

Oi 

1 

T«-IV 

20 

0,62.'> 

SS'ol'' 

■  1,.50 

17*4.5' 

0 , 658 

o 

TR-IV 

20 

0...575 

89'0.5' 

1,50 

l.S*45' 

0,790 

■ 

TR-IV 

20 

0,.7!)0 

87551' 

1  .50 

17*42' 

0,970 

4 

TF-aV 

lo.  20.  30 

0,62.5 

88*51' 

30  MM 

1S*.52' 

0,900 

0 

TR-2V' 

20 

0,.570 

89*05' 

1,50 

0,910 

6 

TB-IV  —  1 

20 

0,600 

90* 

1,50 

T!)"00' 

0,810 

7 

TB-IV  —  ' 

20 

0,62.5 

90* 

1,.‘0 

1!)°26' 

0,8,50 

■S 

T«-:v  ~ 

20 

0,.57,5 

89*05' 

1,50 

0,800 

.  9 

20 

0,600 

90* 

i  ,.50 

19*00' 

0,810 

10 

tbf-Iv  —If' 

20 

0,62.5 

90* 

1,50 

19*20' 

0,8.50 

11 

TR-IV  —  .2‘' 

20 

0,62.5 

90* 

1,50. 

18*50' 

0,920 

12 

TR-IVS  —  2-' 

20 

0,625 

88*,3-t' 

1,50 

US' 25' 

0,827 

lo 

TR-IVS  ~ 

20 

0,625 

.S,S*31' 

1,50 

18*25' 

■  0,910 

:14 

TR-lVS  —  S’-' 

20 

0,625 

88*34' 

1..50 

■  18*25' 

■  1,00 

The  outlet  section  of  the  airfoil  back  is  bent* 
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hj 

Fig.  2.  Pressure  distribution  along  the  airfoil 
In  a  cascade:  a)  for  cascade  No.  2  as  a  function 
of  the  Mach  number  -=  3-7°);  b)  for  cascade  No. 

4  as  a  f'unctlon  of  the  flow-inlet  angle. 
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Fig.  4.  The  spectra  of  flow  In  cascades:  a;  ca 
cade  No.  8,  M  =  1.68,  -  22°30* ;  b)  cascade  No 

9,  HI  ===  1.69,  =  22°30»;  c)  cascade  No.  12,  M  = 


lit  0 


<i) 


Pig.  5.  Iioss  distribution  along  the  height  of  the  cascade 
1)  cascade  No.  4,  b  =  15  mk,  f3^^'==  21°;  b}  cascade  No.  9, 

3  =  20  nmi,  =  20°;  c)  cascade  No.  11,  b  -  20  nun,  f3^  -  20 

i)  cascade  No.  Ip,  b  =  20  mm,  [5  20°. 


aj 


O  3  6  9  13  HM 


Pig.  6.  Distribution  of  flow-outlet  angles  along  the 
height  of  the  cascade:  a)  cascade  No.  4,  b  =  15°  mm, 
=  21° j  b)  cascade  No.  12,  b  =  20  mm,  0^  =  20°. 
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AN  APPROXIMATION  METHOD  FOR  CALCULATING  STEADY-STATE  UNIFORM 
GAS  FLOW  IN  A  TUBE  WITH  EQUILIBRIUM  DISSOCIATION  AND 
IONIZATION,  TAKING  INTO  ACCOUNT  RESISTANCE 
AND  HEAT  TRANSFER 

I.  I.  Suksov 

At  high  temperatures  we  must  take  Into  account  gas  dissociation 
and  Ionization,  which  greatly  complicates  the  calculation  of  gas  flow. 
For  equlllhrium  processes  the  calculations  are  simplified  considerably 
If  we  have  tables  of  the  thermodynamic  functions.  Such  tables  are 
available  for  air  [l,  2].  The  handbook  "Physical  Gasdynamics"  [5] 
contains.  In  addition  to  the  theoretical  bases  for  the  compilation  of 
these  tables,  the  solution  of  certain  gasdynamlc  problems  for  dissoci¬ 
ated  and  Ionized  air.  Individual  problems  of  gasdynamics,  taking 
Into  account  the  Influence  of  dissociation  and  ionization  of  the  air, 
have  been  treated  previously  [4]. 

On  the  basis  of  the  tables  [l,  2]  data  have  been  obtained  from  a 
calculation  of  one -dimensional  Isentroplc  air  flow  with  equlllbrlurt 
dissociation  and  Ionization  [12]. 

In  this  work  we  give  an  approximate  graphoanalytlc  method  which 
makes  it  possible  to  calculate  one -dimensional  Isentroplc  flow  In 
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tubes,  l.e.,  to  take  Into  account  the  influence  of  resistance  and  heat 
transfer. 

Yu.  I.  Balakleyevskiy  helped  in  the  calculations  used  to  develop 
this  method. 


1.  An  Approximate  Solution  to  the  Equations 


Let  us  write  the  equations  of  steady-state  one -dimensional  move¬ 
ment  of  a  viscous  compressible  gas  in  a  tube  during  constant  flow  with 
removal  of  heat  from  the  gas  [5]: 


flow  equation: 


energy  equation: 


Fgow  =  G  =  const. 


/gdh  -i-wcIw-T  Jgdqa  —  0, 


equation  of  the  first  law  of  thermodynamics; 


(1.1) 


(1.2) 


(1.3) 


where  E  is  the  cross-sectional  area  of  the  tube  (m'^);  £  is  the  accelera- 
tlon  of  gravity  (m/sec  );  p  is  the  density  (kg-sec  /m^);  w  is  the 
velocity  (m/sec )j  G  Is  the  weight  flow  (kg/sec);  I  =  427  kg-m/kcal  is 
the  mechanical  equivalent  of  heat;  h  is  the  enthalpy  (kcal/kg);  dq  Is 
the  elementary  amount  of  heat  (kcal/kg);  dq..  Is  the  elementary  amount 

a. 

of  heat  released  from  1  kg  of  gas  (kcal/kg);  £  is  pressure  (kg/m  ); 

2  2 

and  dZp  is  the  elementary  work  of  drag  (ra  /sec  ). 


On  the  basis  of  (I.5)  we  get  an  equation  for  entropy  £: 


The  value  of  di^  is  defined  by  the  expression 

ai’  dx 


(1.4) 


(1.5) 
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where  ^  Is  the  resistance  coefficient;  x  Is  the  coordinate  read  along 
the  tube  axis  (m);  D  Is  the  tube  diameter  (m),;  for  tubes  of  noncircular 
cross  section^  D  Is  understood  to  be  the  hydraulic  diameter. 

The  elementary  heat 


where 


='7*-r';rad. 


(1.6) 

(1.7) 


O 

dPg  Is  an  area  element  of  the  side  of  the  tube  (m  );  Is  heat  flow 

o  2 

to  the  wall  (kcal/m  sec);  is  convective  heat  flow  (kcal/m  sec);  and 

2 

'^rad  radiation  heat  flow  (kcal/m  sec). 

The  following  relations  hold: 


(1.8) 

Cjh  =  «/,  (/^i 

(1.9) 

(1.10) 

, ^-0,001361 

(1.11) 

2 

where  Is  the  local  heat-transfer  coefficient  (kg/m  sec);  = 

=  0.  5(1  +  ^h  ^  local  Stanton  number; 

^wall  degree  of  blackness  of  the  wall;  e  Is  the  degree  of  blacl 

ness  of  the  gas  at  temperature  T;  A  Is  the  absorptivity  of  the  gas  at 


temperature  T  (assimed  equal  to  the  value  of  e  at  temperature  T  ). 

w  w 


The  subscripts  1  and  w  Indicate  conditions  of  a  heat-insulated 
and  cooled  wall,  respectively. 


3  Ions 


Taking  Into  account  (i.5)-(l.ll)  we  get  for  dq„  and  ds  the  expres- 

GL 

(l.lj) 

The  enthalpy  h^^^  Is  defined  by  the  formula 
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(1.14) 


•  A  +  r^,,- 


2/^ 


where  r^  is  the  recovery  factor. 

Let  us  examine  the  problem  of  determining  the  velocity  w,  tempera¬ 
ture  T,  pressure  £,  and  density  p  of  a  hlgh-temperature  gas  in  cross 
sections  of  a  tube  In  presence  of  dissociation  and  Ionization,  consider¬ 
ing  these  processes  to  be  equilibrium  processes.  The  starting  data 
are:  gas  flow  G,  temperature  Tqq  and  pressure  p^^  of  Isentroplc  stag¬ 
nation  in  cross  section  x  =  Xq,  and  also  the  geometry  of  the  tube.  In 
addition,  we  should  know  the  dependences  of  the  dimensionless  values 
rj^,  and  on  the  similarity  parameters;  the  degree  of  blackness  of 
the  gas  e;  the  degree  of  blackness  of  the  surface  of  the  tube  walls 
^wall^  and  we  should  have  the  distribution  of  enthalpy  or  temperature 
on  the  wall  for  the  length  of  the  tube. 

As  will  be  shown,  with  tables  of  the  thennodynamlc  functions 
similar  to  those  for  air  [1,  2],  the  problem  can  be  solved  relatively 
simply  by  means  of  graphoanalysis. 

Using  Eqs.  (l.l),  (1.2),  (1.4),  and  (l.5),  the  system  of  equations 
for  the  given  problem  can  be  written: 


w  =  V2lg  / A,„  -  -  A  -  -  A/;, 

*  ~  *00 
G  —  '(wF, 


(1.15) 

(1.16) 

(1.17) 

(1.18) 


where  y  =  gp;  a  Is  the  speed  of  sound; 


iix_\ 

2ig  ■  o  y  ■ 


The  value  of  dq^  Is  determined  from  (l. 12).  The  subscript 

a 

Indicates,  as  before,  stagnation  parameters  when  x  =  Xq. 

Let  us  first  examine  flow  without  resistance  and  heat 
(dlj,  =  dq^  O),  i.e.,  Isentroplc  flow. 


(1.19) 

"00" 

transfer 
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For  this  case,  Eqs.  (l. 15)-(1. l8)  assume  the  form 


=  /2/^  /A,,,,  ’  A, 

(1.20) 

s  ==s„„  =  const, 

(1.21) 

T 

(1.22) 

i 

■w  =  aM. 

(1.23) 

Knowing  Tq^^' and  p^Q  let  us  determine,  from  the  tables  of  thermodynamic 
functions,  the  enthalpy  hQQ  and  the  entropy  Sqq.  For  a  number  of  values 
of  pressure  p  <  Pqq  and  when  s  =  Sqq,  we  determine,  from  the  tables, 
the  corresponding  values  of  temperature  T,  enthalpy  h,  speed  of  sound 
a,  and  the  value  of  y.  Using  Eq.  (l.20)  let  us  construct  a  graph  of 
w  vs.  h;  using  (I.23)  let  us  construct  w  vs.  h  for  a  ramiber  of  given- 
values  of  M  (a  set  of  curves).  The  Intersection  of  the  curves  deter- 
mines  the  values  of  w  and  h  corresponding  to  the  selected  M  values. 

Then  (graphically  or  by  Interpolation  using  the  tables)  we  determine 
the  desired  values  of  T,  £,  and  p,  and  from  Formula  (1.22)  we  calculate 
the  values  of  the  cross-sectional  area  F  (if  £  Is  given,  we  determine 
the  flow  from  (1.22)). 

The  solution  for  Isentroplc  flow  Is  the  first  approximation  when 
solving  the  problem  of  gas  flow  In  a  tube. 

This  method  Just  examined  differs  somewhat  from  that  used  previously 
[12].  We  note  that  the  graphic  data  [12]  In  many  cases  are  sufficient 
to  determine  all  the  parameters  of  Isentroplc  air  flow  and  for  a  logi¬ 
cal  division  of  the  tube  Into  parts  In  order  to  take  Into  account  the  ' 
Influence  of  resistance  and  heat  transfer. 

Now  let  us  examine  the  solution  of  Eqs.  (l. 15)-(l. I8)  for  the 
general  case.  Let  us  divide  the  tube  into  parts  and  examine  the  k-th 
section.  The  values  at  the  beginning  and  end  of  this  section  will  be 
designated  by  the  subscripts  (k-l)  and  k,  respectively. 
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On  the  basis  of  Eqs.  (l. 15)-(i- l8)  we  have 


2lg  V Ao.  —  A*  —  AA*. 

s*  =  s*-i 


W.  =- 


IkFk 

«k  \ 


(1.24) 

(1.25) 

(1.26) 

(1.27) 


where  hQ  =  w^_^/2Ig  is  stagnation  enthalpy  in  the  (5c-l)-th 

section, 

■'  b  "*4 

As*  =  J  -{dq,  -  --  —  j  .  (  i.  28 ) 

•>>-1  -^*-1 

where  dq  is  taken  on  the  basis  of  (1.12). 

3. 

The  Integrals  in  each  section  are  defined  approximately  by  the 
system 

/*==J  (1.29) 

jr*_l 

where  the  subscript  "am"  Indicates  the  taking  of  the  arithmetic  mean. 

Ax  =  jc*-x*_,. 

Using  the  known  parameters  In  the  (k-l)-th  section,  the  gas  para¬ 
meters  In  the  k-th  section  are  determined  as  follows. 

First  we  find  the  gas  parameters  In  the  k-th  section,  considering 
that  the  flow  In  the  examined  section  Is  Isentroplc  (these  parameters 
are  given  the  subscript  "kl").  Knowing  the  parameters  with  the  sub¬ 
scripts  "k-l"  and  "kl"  let  us  calculate  the  Integrals  according  to 
system  (1.29);  let  us  determine  the  entropy  according  to  (1.25). 
Assigning,  when  s  =  Sj^,  several  test  values  to  Pj^  close  to  Pj^^  and 
determining  hj^  and  by  Interpolation  of  the  tables,  we  calculate  the 
corresponding  values  of  Wj^  from  Formulas  (1.24)  and  (1.26)  (the  values 
of  Fj^  are  the  same  as  those  obtained  earlier  for  Isentroplc  flow).  Let 
us  construct  the  curves  Wj^  =  i*^(h^)  and  Intersection 
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of  these  curves  gives  the  desired  values  of  and  hj^.  Then  we  deter¬ 
mine  the  parameters  Tj^,  a^^,  and  Pj^,  and  calculate  the  number  M 
from  (1.27).  Thus  is  calculated  flow  in  a  tube^  beginning  with  the 
first  section.. 

In  this  calculation  method,  as  the  thermodynamic  magnitudes  which 
uniquely  define  all  parameters  of  gas  flow  in  atuhe,  we  selected 
enthalpy  h  and  entropy  Let  us  note  that  it  becomes  ambiguous  to 
determine  the  gas  parameters  using  total  enthalpy  hQ  ane  entropy  s  = 

=  Sq  in  the  case  of  Isentroplc  flow. 

2.  Data  for  Approximate  Evaluations 

When  using  this  approximation  calculation  method  we  must  know, 
as  has  been  said,  the  dependences  of  the  values  rj^,  and  on  the 
similarity  parameters,  and  the  dependence  of  the  degree  of  blackness 
of  the  gas  e  on  temperature,  pressure,  and  path  length  for  the  dis¬ 
sociated  gas.  As  far  as  we  know,  there  are  no  such  data  for  tubes. 

For  effusers  and  tubes  of  constant  cross  section,  in  particular  for 
nozzles  and  the  working  part  of  thermal  wind  tunnels,  we  can  consider 
that  the  resistance  is  exhausted  by  friction,  l.e., 

^  =  49.  (2.1) 

where  c„  =  2t  /pw  is  the  local  friction  factor,  and  t  Is  frictional 

X  w  w 

stress  on  the  wall  (kg/ra  ). 

In  these  cases,  if  the  botmdary  layer  occupies  a  small  portion  of 
the  cross  section,  the  use  of  rules  of  friction  and  heat  transfer  for 
a  flat  plate  in  the  calculations  will  give  a  good  approximation.  How¬ 
ever,  even  for  a  flat  plate  we  do  not  have  available  systematized  data 
on  friction  and  heat  transfer  for  a  sufficiently  broad  range  of  changes 
in  the  similarity  criteria. 


To  obtain  approximate  evaluations  of  the  Influence  of  convection 
heat  transfer  and  friction  we  can^  until  more  reliable  data  are  forth¬ 
coming,  use  formulas  In  which  account  Is  taken  of  the  variability  of 
the  physical  parameters  across  the  boundary  layer  and  the  basic  effect 
of  dissociation  [6]. 

For  a  laminar  boundary  layer 


where 


9. ■  R<-'  =  0.(>6.! (^)\  (-^y^ , 

T 

C),,  Re"/  =  0,332  Pr,,  )'=  , 


_ .sIm'  c  _  _ _ 1 


(2.2) 


(2.5) 


h^  Is  determined  from  Formula  (l.l4)  where  r  =  0.84  and  Pr^  =  0.7I  = 


=  const.  For  a  turbulent  boundary  layer 

Cy-a  Re^;’  =  0,058  •  (1  -|-O,2rA/'0"’", 

C,„„Re«/  -=  0,029  P,-/"  ■  ( H'  0,2rM')"’‘ ' , 


(2.4) 

(2.5) 


where  r  =  O.89.  The  Itinctlon  h(T)  is  given  In  Fig.  1. 

The  subscript  Indicates  parameters  pertaining  to  maximum 


enthalpy  In  boundary  layer  h,  whereupon  we  have  the  relation 

A,  — fta 


ha  —  ha 


=  0,25 


=  --HL 

2/i'-  ' 


(2.6) 


Exceptions  to  this  formula:  If  co  <  1  -  h^h,  h;,^  =  h;  If  h^h^  > 

=  V 

The  formulas  for  the  transition  from  c^  and  to  c^  and  C^: 


Pi- 


?  P 


Figure  2  gives  the  coefficient  of  viscosity  p,  vs.  h  for  T  <  i875°K. 
When  T  >  200°K  we  can  use  the  familiar  graphic  function  m.(T,  P)  [11]. 
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Let  us  note  that  the  last  term  In  (2.2)  and  (2.5)  can  he  dis¬ 
regarded  because  of  the  small  e^qsonent. 

¥e  should  bear  In  mind  that  the  anproxlmate  nature  of  data  on 
convection  heat  transfer  and  friction  Introduce  no  noticeable  error 
Into  the  flow  parameters,  since  the  estimated  Influence  Is  relatively 
slight,  even  In  a  turbulent  regime. 

The  values  of  ^^®cessary  for  calculating  radiation  heat  trans 

fer  are  familiar  for  a  great  many  materials  [7,  81.  The  accessible 
literature  contains  no  data  on  the  degree  of  blackness  of  a  gas  e  at 
temperatures  higher  than  2000°K.  For  temperatures  T  <  2000*^K  we  have 
data  on  gas  radiation  [7l»  In  which  Is  noted  the  negligible  radiation 
absorptivity  of  monatomic  and  diatomic  gases.  We  also  know  [9]  that 
at  low  pressures  (p  <  0.1-  atm)  the  role  of  atmospheric  radiation  In 
the  general  heat  balance  Is  slight,  even  at  very  high  temperatures 
(8000-10, 000°K).  Prom  physical  concepts  It  follows  that  with  increas¬ 
ing  pressure  and  Increasing  beam  path  length  the  radiation  absorptivity 
of  a  gas  Increases  [7],  l.e..  Its  degree  of  blackness  Increases,  tend¬ 
ing  toward  a  limiting  value  e  =  1  for  an  absolutely  black  body.  Such 
conditions  can  be  realized,  e.g.,  with  high-voltage  arc  discharge; 
then  It  Is  assumed  that  the  arc  radiates  as  an  absolutely  black  body 
[10]. 

5.  Calculation  Results 

/ 

To  extlmate  the  combines  Influence  of  convection  heat  transfer, 
radiation,  and  friction  we  calculated  the  air  friction  In  a  section 
of  a  round  cylindrical  tube  with  diameter  d  ■  20  mm  and  length  1-50 
mm  (l/d  =  2.5)*  The  flow  regime  was  assumed  to  be  turbulent.  We  used 
the  following  data:  stagnation  temperature  and  pressure  In  the  Initial 
section  Tqq  =  12,000°K,  Pqq  =  200  atm;  M  =  1;  cooled-wall  temperature 
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=  9Q0°K;  the  value  0.5(1  +  degree  of  hlackness  of  the 

air  e  w  1.  The  section  was  not  divided  into  parts,  since  the  tube  was 
cylindrical  and  the  ratio  l/d  was  small.  After  determining  the  para¬ 
meters  of  isentroplc  flow  we  calculated  the  drop  In  total  enthalpy 
Ah  and  entropy  As  from  the  formulas 


A//  =  AA*  +  l/!i,ad.  -^5  = 


A/t*  =  0,116  Pr^-o."  ^  0,2rM=f"  (h^  \  -  hj -j 


AA. 


"r«i=  0.0054+  ( -V-'-  L 
\100/  7a.  (I 

AA,  =  0,1 10  +  0,2rAl-)''' 


II  l 
Ig  d 


The  derived  working  formulas  or  Ah,  ,  Ah  ,,  and  Ah-  were  obtained 

k  rad  f 

on  the  basis  of  relations  (1.12),  (I.19),  (2.1),  (2.4),  (2.5),  and  (2.7). 
When  calculating  Ah^^^  the  terra  A(T^IOO)^  was  discarded. 

We  obtained  the  following  values:  Ahj^  =  514,  Ah^^^  =■  390,  Ah^  = 

=  67,  Ah  =  904,  As  =  0.0758. 

Then,  using  the  given  method  we  determined  the  flow  parameters 
at  the  end  of  the  tube. 

Figure  3  and  4  gives  the  graphic  determination  of  the  parameters 
w,  h,  T,  p,  and  p  for  Isentroplc  flow  and  at  the  end  of  the  tube, 
respectively  (the  velocities  w^  were  calculated  from  (l.20)  and  (1.24)). 
The  table  gives  the  calculation  results. 

This  example  of  a  calculation  shows  that  heat  transfer  and  fric¬ 
tion  can  have  a  substantial  Influence  on  the  parameters  of  the  flow  of 
a  high- temperature  gas  In  a  tube.  This  can  happen,  e.g..  In  the  Initial 
sections  of  nozzles  of  supersonic  hyperthermal  tubes. 
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